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aeMt^tT«TiSaervKJich  has  been  developed  to  pre- 
dict the  thermal  response  of  deep  sea  diverb  working  at  depths  as  great  as 
1000  feet.  Under  these  conditions,  divers  are  subjected  to  extreme  thermal 
stress  because  the  water  temperature  is  approximately  0°C;  the  gas  used  for 
breathing  and  in  the  suit  is  helium  which  has  a high  thermal  conductivity. 
Increased  mass  flow  rate  through  the  respiratory  tract  promotes  heat  loss  when 
the  breathing  gas  is  not  warmed.  Since  the  experimental- evaluation  of  new 
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(continued) 

systems  designed  to  protect  divers  is  a difficult,  time  consuming  process 
it  is  worthwhile  to  have  a reliable  mathematical  model  for  predicting  diver 
performance  under  various  conditions.  8 

...  "fjhematical  model  described  in  this  paper  subdivides  the  human 

into  i5  cyUndricai  elements  representing  the  head,  thorax,  abdomen,  and 
proximal,  medial,  and  distal  segments  of  each  arm  and  leg.  Within  each 
element  the  transient-state  heat  conduction  equation  is  solved  to  obtain 
temperature  profiles.  Mass  balances  are  also  computed  for  oxygen,  carbon 
oxide  and  lactic  acid.  The  program  permits  one  to  specify  various  kinds 
of  wet  and  dry  suits  with  either  open  or  closed  circuit  heating. 


Representative  results  are  presented. 
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IN  TROD"  Us 

On  June  5-6,  IS  5 a ■ aa-trcial  diving  te  iu>  descended  to 
a depth  < i 30. » on  ters  : tne  Lay  ot  Laburador  a.  <i  helped  to 
re-e^iaMlsh  a drilling  . erat Ion  that  had  been  lnterupted 
b • a violent  storm.  While  this  was  not  a tvoicn.  . ummt  r-ial 
d Ivin*  job.  f wa-  indicative  ot  things  to  come-  Today, 
dives  to  depths  ot  /.GO  meters  are  becoming  commonplace  in 
the  Nortl  Sea.  and  'cadiug  comaer  i»I  diving  firms  are  devel- 
oping systems  ' i working  at  greater  depths  as  drilling 
operation*.  mo'v  further  offshore. 

Hear  loss  t (-presents  i serious  problem  for  d Ivors 
under  then*,  condit  • « u» . Th*-  severe  thermal  stress  Imposed 
on  divers  n i maw  r ft  Ion  in  cold  water  is  compounded  by  factors 
relate^  ch.  hvperbaric  environment.  At  depths  greater 
than  : 0O  ruetfc.  the,  breathing  gas  is  usuall;  a mixture  of 
helium  and  xygen,  w!.i'h  ha*  a thermal  conductivity  approx- 
imacel  - ..ik  tinwrs  that  or  air  This  causes  a marked  de- 
crease • r .j]  insulation  of  such  garments  as  woolen 

undone  u me!  ! aaed  neoprene  outer  suits.  Furthermore, 
cl-.-sed  - ' f..med  T.uteilals  tend  to  compress  and  lose  their 

insular!  up.  .-.hi  lit-,  it  depth.  Heat  lots  through  the  respi- 
ratory t.a  t 'l.o  Increases  iwing  to  an  Increase  in  gas 
density  wl-  ii'  •- easing  depth.  'f  the  breatr.ing  gas  is  not 
heated,  heat  through  the  respiratory  tract  will  be 

approximate  .<  t.,  the  rate  f metabolic  heat  generation 

at  a de: * • '00  m*  rets,  and  since  the  ventilation  rate  is 

roughlv  r.  , nal  r the  m.  tabolic  rate  for  moderate  ex- 
ercise levels  i will  be  true  even  when  the  diver  is 

working  < : ) 

Va-lous  sy  .r-  , have  been  developed  to  provide  thermal 

prote  tion  f r ..  or,.  At  moderate  depths  ro-nmercial  divers 
often  wear  a - ? - i t flooded  wifi,  warm  water  supplied  from 
the  surface  Military  divers  engaged  in  relatively  shallow 
clandestine  missions  must  rely  on  self-contained  systems, 
usually  consisting  of  only  a wet  or  dry  suit  without  active 
heating.  Chemical  and  electrical  energy  sources  are  under 
development  for  use  when  passive  suits  provide  inadequate 
protection . 

Experimental  testing  of  such  systems  in  hyperbaric 
chambers  whi^h  approximate  conditions  encountered  in  the 
open  sea  Is  very  expensive  and  time  consuming  Therefore, 
it  is  desirable  to  have  a mathematical  model  that  can  be 
used  both  to  design  new  systems  and  to  Interpret  test  and 
field  measurements.  The  purpose  of  this  paper  is  to  describe 
such  a model  and  discuss  a few  results  obtained  using  the 
model . 

DESCRIPTION  OF  A HEM  AND  MASS  TRANSPORT  MODEL  FOR  THE  HUMAN 


General  Features 

The  phvsical  system  on  which  the  equations  of  change 
are  based  is  shown  in  Fig.  1.  It  consists  of  a number  of 
cylindrical  dements  representing  longitudinal  .segments  of 
the  head,  trunk,  arms,  and  legs.  Each  segment,  consisting 
of  a conglomeration  of  tissue,  hone,  fat,  and  skin,  has  a 
vascular  system  composed  of  three  parts  representing  the 
arteries,  veins,  and  capillaries.  The  circulatory  path  is 
faithfully  reproduced  in  the  sense  that  a tracer  material 
introduced  info  an  artery  lr.  a given  segment  will  flow  both 
into  rhe  capillaries  of  that  segment  and  into  the  arteries 
of  more  distal  segments.  Blood  leaving  the  capillary  bed 
flows  into  a vein  where  it  is  mixed  with  venous  blood  from 
more  distal  segments.  The  mixed  venous  stream  at  the  heart 
flows  into  the  pulmonary  artery;  exchange  of  both  hear  and 
mass  occurs  in  the  lungs. 


Within  a given  segment,  physi  al  paiaae’ers,  such  is 
density  and  specific  heat,  vary  with  radial  position,  wiile 
dependent  variables,  such  as  temperature,  O tension,  meta- 
bolic rate,  and  blood  perfusion  rate,  depend  on  both  position 
and  time.  Certain  variables,  perfusion  rate  and  metabolic 
rate,  for  instance,  are  defined  by  physiological  control 
equations.  Others,  such  as  temperature  and  0.  tension,  are 
dependent  variables  whose  immediate  values  depend  on  tie  pre- 
vious history  of  the  subject.  Using  this  model  one  an  start 
with  specified  initial  conditions  and  follow  changes  that 
occur  in  tesponse  to  various  imposed  conditions. 


Equat lorn,  of  Change  i r Energy 


Development  of  the  beat  transfer  equat  ' = will  be  dis- 
cussed in  some  detail.  Since  the  equal  tens  for  mass  transfer 
are  based  on  analogous  oncepts,  they  will  be  stated  more 
con.  isely. 

Heat  which  is  generated  by  metabolic  reactions  is  either 
stored  in  the  element,  carried  away  bv  circulating  bio. d,  or 
conducted  to  the  surf  a e where  it  is  transferred  to  the  en- 
vironment. the  basic  equation  for  the  thermal  analysis  is 
the  heat  conduction  equation  given  below: 
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- instantaneous  temperature  of  tie  tissue,  bone, 
or  vise  eta  at  a distance  r from  the  axis  of 
the  i-^1  element, 
density  ->f  tissue, 

specific  heat  of  tissue, 

effective  thermal  conductivity  of  tissue, 
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metabolic  heat  generation  rate  per  unit 
volume , 

product  of  the  mass  flow-rate  and  specific 
heat, of  blood  entering  the  capillary  bed* 
per  unit  volume, 

» heat  transfer  coefficient  between  arteries 
and  tissue  per  unit  volume, 

• heat  transfer  coefficient  between  veins  and 
tissue  per  unit  volume, 
temperature  of  arterial  blood, 

temperature  of  venous  blood. 

It  should  he  observed  that  this  form  of  the  heat  conduction 
equation  is  applicable  only  to  an  axially  symmetrical  system 
In  which  longitudinal  conduction  of  heat  is  negligible.  If 
the  subject  is  moving  so  that  there  Is  a uniform  flow  of  air 
or  water  around  each  of  the  elements,  the  analysis  should 
apply.  H.  H.  Pennes  (2)  has  shown  that  longitudinal  con- 
duction in  the  arms  is  relatively  unimportant,  which  should 
be  true  also  In  the  legs  (3),  but  not  in  the  head. 


T ft) 


Source  terms  on  the  right-hand  side  of  Eq . 1 represent 
the  rate  of  heat  generation  bv  metabolic  reactions,  the  net 
rate  at  which  heat  is  transferred  from  artetial  blood  to 
tissue,  and  the  rate  at  which  heat  is  translerred  from  ven- 
ous blood  to  tissue.  It  is  assumed  that  there  is  perfect 
heat  transfer  between  blood  in  the  capillaries  and  neighboring 
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tissue;  i.e.,  the  temperature  of  blood  leaving  a apiliary 
bed  is  equal  to  the  temperature  of  the  neighboring  tissue. 

As  a first  approxlmat ion , It  is  assumed  also  that  the  rate  of 
heat  transfer  ft  os  blood  in  large  vessels  to  neighboiing 
tissue  is  proportional  tj  the  difference  between  the  blood 
and  tissue  temperatures  The  proportionality  factois  are 

called  H . for  arteries  and  H , for  veins, 
al  v i 

Since  the  temperature  of  blood  in  large  vessels  changes 
with  time,  it  Is  necessary  to  write  two  more  thermal  energy 
balances  In  formulating  an  equation  for  arterial  blood,  we 
assume  that  arteries  in  the  i*Jl  element  torn  a pool  having  a 
uniform  r eap-.rature , T . The  rate  of  accumulation  of  ther- 
mal energy  in  this  reservoir  la  equal  to  the  sum  of  the  net 
rate  at  which  heat  is  arned  into  the  pool  by  flowing  blood, 
the  rate  at  wnlth  heat  Is  transferred  iron  neighboring 
tissue  to  blood  in  the  pool,  and  the  rate  at  which  heac  is 
transferred  directly  from  the  venous  pool  to  the  arterial 
pool  owing  to  the  proximity  of  certain  arteries  and  veins. 
This  equality  Is  expressed  mathematically  by  the  following 
equation . 


in  which 


product  of  the  mass  flow  rate  and  specific 
heat  for  arterial  blood  flowing  into 
capillaries , 

product  of  the  mass  flow  rale  and  specific 
heat  for  venous  blood  flowing  into  pulmonary 
capillaries . 


Equations  2 and  3 must  be  modified  to  take  cognizance  of  the 
fact  that  venous  blood  flows  into  the  pulmonary  capillaries 
which  in  turn  empty  into  the  arterial  pool: 


(t)  - temperature  of  tne  blood  entering  the  arterial 
pool , 

■ mass  of  the  blood  contained  in  the  arteiial  pool 
ot  Che  iii?  element, 
specific  heat  of  blood. 

It)  - product  of  the  mast,  tlow  rate  and  specific 
he..t  for  blood  entering  the  arterial  pool, 

• length  of  the  t^l)  clement, 

t “ hear  transfer  coefficient  for  direct  transfer 
between  large  arteries  and  veins. 


rate  at  which  heat  is  tiansferred  from  venous 
blood  in  the  thorax  to  air  in  the  respiratory 
tract, 

rate  at  which  venous  blood  flows  from  the 
element  Into  the  venous  pool  in  the  thorax  » 

0 j(t)  for  those  elements  which  are  connected 
J to  the  thoracic  segment. 


The  total  rate  of  heat  loss  through  the  respiratory  tract 
depends  on  the  ventilation  rate  and  the  temperature  and  humid- 
ity of  inspired  air.  In  this  analysis,  it  is  assumed  that 
10  percent  of  tne  heat  loss  through  the  respiratory  tract 
comes  from  the  arterial  pool  in  the  head,  10  percent  from  the 
venous  pool  in  the  head,  and  80  percent  from  the  venous  pool 
in  the  thorax. 


The  corresponding  equation  for  the  venous  pool  is 


Before  these  equations  can  be  solved  uniquely,  certain 
constraining  conditions  must  be  specified.  Some  of  these  take 
the  form  of  Initial  conditions  which  specify  all  temperatures 
at  the  instant  the  transient  begins: 


product  of  mass  and  specific  heat  for  venous 
blood  flowing  into  the  venous  pool  of  the 
iil*  element  from  the  (A--  element. 


Since  It  is  assumed  throughout  this  analysis  that  and 
M are  constants, 


Also  needed  are  boundary  conditions  which  relate  the  subject 
to  his  environment.  These  take  the  form: 


The  aquation  for  the  venous  temperature  in  the  abdominal 
section  is  slightly  different  from  F.q . 3 bet  ause  a vein 
from  each  leg  flows  into  this  section.  It  is  also  necessary 
to  modify  the  equarlons  for  the  thoracic  section  because  all 
venous  streams  terminate  and  arterial  streams  originate  in 
this  section.  It  is  assumed  that  the  temperature  of  blood 
entering  the  pulmonary  capillaries  is  equal  to  the  "cup 
mixing"  mean  temperature  of  venous  streams  entering  the  right 
ventricle.  This  necessitates  a change  in  Eq . 1 because  the 
temperature  of  venous  blood  entering  the  pulmonary  capil- 
laries is  different  from  the  temperature  of  arterial  blood 
entering  more  superficial  capillaries  of  the  thorax: 


1 effective  environmental  temperature 
rate  of  heat  loss  by  evaporation. 


in  which 


heat  transfer  coefficient  for  convection 
heat  transfer  coefficient  for  radiation. 


For  unclothed  regions  of  the  body,  the  surface  is  skin 
but  for  clothed  regions,  it  is  the  outer  surface  of  the 


garment.  Since  physical  proper  tie;-  are  permitted  to  vary 
with  radial  position  in  this  model,  garments  of  considerable 
complexity  can  be  analyzed.  One  can  include  a layer  of  g<»s 
or  water  between  skin  and  the  garment  by  assigning  suitable 
values  for  the  thermal  conductivity,  density,  and  specific 
heat.  It  is  also  possible  to  include  an  unde « garment  which 
contains  tubing  through  which  warned  or  cooled  liquic  passes 


Equations  ot  Change  t o r Mas s 


Corresponding  equations  for  mass  transfer  are  similar 
to  those  presented  above  For  accumulation  of  particular 
component  in  tissue  and  blood,  we  have  the  following  equa- 
tions. 

oC 


(13) 


and 


?V  5, 

b 9t  11  3x 


(14) 


in  which 


(0.24  £ +•  C > UN>  ,i 

<T  dt  1 q <C02.i  ‘ C02.i-1>  * 4~*t 


If  0 ' CQ  t « .046 


(1?) 


e dC0,i  t 

r— fr  + q <co2.i-co2.i-i>-r  h 


If  .046  C, 


0 , i . 


The  subscript  b has  been  dropped,  but  it  is  to  be  understood 
in  the  remainder  of  the  p£ger  that  Cy  refers  to  the  oxygen 
content  of  blood  in  the  i—  segment  2’  of  the  capillary 
(1  < i 4) . For  the  first  segment,  ^ j refers  to  the 

oxygen  content  of  arterial  blood  2*  entering  the 

capi 1 larv. 


e,  and  e b • volumetric  fractions  of  tissue  and  blood, 
respectively, 

and  C.  - concentration  in  tissue  and  blood,  re- 
C D , 

spectively, 

q » flux  of  blood  through  the  capillaries  ir.  a unit 
area , 

x “ distance  measured  along  a capillary, 

l'tb  * rate  of  exchange  between  tissue  and  capillary 
per  unit  volume, 

R • rate  of  production  per  unit  volume. 

Although  trie  basic  equations  for  all  components  have  the 
same  form,  the  approxi mac  ions  which  are  appropriate  in 
dealing  with  the  distribution  of  material  between  tissue  and 
blood  depend  on  the  particular  omponent  involved. 

When  dealing  with  oxygen,  it  is  assumed  that  equili- 
brium exists  between  oxygen  in  myoglobin  and  hemoglobin. 
Typical  functional  relationships  between  oxygen  content  and 
oxygen  tension  for  myoglobin  and  hemoglobin  are  shown  in 
Pig.  2.  Since  myoglobin  is  essentially  saturated  at  oxygen 
tensions  greater  than  20  Torr,  the  amount  of  oxygen  stored 
in  myoglobin  remains  constant  under  notmal  conditions.  These 
stores  arc  released  when  the  oxygen  tension  falls  below  20 
Torr  The  dissociation  curve  lor  myoglobin  can  be  approxi- 
mated by  the  dashed  curve  without  introducing  appreciable 
error.  Then  equilibrium  between  oxygen  in  myoglobin  and 
hemoglobin  implies  that 

Cn  * 0.24  C , if  0 <-  C .046  volume  fraction 

02*  ~ °2b  _ 

- 0.011  if  .046  < . . (15) 

°2b 


Whenever  the  max. mum  rate  of  ox/ger.  transport  to  muscle 
exceeds  the  demand,  the  rate  of  utilization  is  determined  by 
the  demand.  During  vigorous  exercise,  however,  demand  ex- 
ceeds supply  and  anaerobic  reactions  occur,  if  falls 
below  a critical  value  in  any  capillary  segment,  2 it  is 
assumed  that  the  oxygen  supply  is  Inadequate  and  anaerobic 
glycolysis  occurs  with  the  production  of  la-  tic  acid.  The 
critical  oxygen  concentration  was  arbitrarily  set  at  a con- 
stant value  of  7.5  volume  percent  although  it  should  depend  on 
the  density  of  capillaries  and  the  local  metabolic  rate,  both 
of  which  change  during  exercise. 

Since  the  set  of  rea<  t ions  involved  in  metabolism  is 
very  romp  lex.  w*.  found  it  necessary  to  adopt  a greatly  simpli- 
fied metabolic  scheme  whiih  Joes,  nevertheless,  include  sev- 
eral of  the  more  important  factors.  Aerobic  synthesis  of 
ATP  using  glucose  as  th«  substrate  proceeds  according  to  the 
following  overall  reaction. 

C,H.,0,  +60,+  38  ADP  ♦ 18  P 

6 12  6 2 i 

(glucose) 

* 6 CO 2 + 6 H2O  + 38  ATP  (18) 

In  the  absence  of  oxygen,  synthesis  of  ATP  proceeds  as  f < . 1 — 

1 ows . 

C.H . _0  e 2 ADP  + 2 P . 
o IJ  n 1 

(glucose) 

; 2 C,H,0,  + 2 AiP  (19) 

i 6 J 

(lactic  acid) 


When  Eqs.  13  and  14  are  added,  one  obtains 


dco2b  *0,b 

(0.24  r 4-  e.  ) -3- + q — i-  - CM 


3Co.b  aco2b 

1 9 — • e R If  .046  *■  C. 


tb  3t 


(16) 


In  this  case,  R is  a negative  quantity  which  accounts  for 
depletion  of  oxygen  by  metabolic  reactions. 

Allowance  is  made  for  the  reduction  in  oxygen  tension 
along  a capillary  by  subdividing  it  into  four  section*,  each 
treated  as  a well-mixed  volume.  Then  Eq.  16  assumes  the 
form 


Lactic  acid  produced  in  muscle  is  transferred  to  blood  in  the 
capillaries  and  distributed  throughout  the  body.  It  can  be 
converted  to  glycogen  in  the  liver  according  to  an  aerobic 
mechanism  which  oxidizes  one  mole  of  lactic  acid  to  carbon 
dioxide  and  water  for  every  four  moles  converted  Into 
glycogen,  or  it  can  be  cataholized  in  muscle. 

We  assume  that  glucose  is  always  available  at  the  rate 
required.  Hence,  it  is  onlv  necessary  to  account  for  three 
components  - oxvgen,  carbon  dioxide,  and  lactic  acid. 
Specification  of  the  rate  of  energy  release  assuming  that 
only  aerobic  processes  are  Involved  defines  the  demand  for 
oxygen.  If  the  rate  of  oxygen  supply  as  defined  above  is 
not  adequate  to  meet  the  demand,  the  deficit  will  be  pro- 
vided by  anaerobic  reactions.  However,  19  moles  of  glucose 
are  consumed  In  the  anaerobic  production  of  the  same  amount 
of  ATP  that  is  produced  using  one  mole  of  glucose  in  the 
aerobic  process.  Knowing  this  one  can  compute  production 
rates  lor  carbon  dioxide,  lactic  acid,  and  heat  for  given 
values  of  demand  and  supply. 

Equations  of  change  for  carbon  dioxide  are  derived 
using  the  same  approach.  It  is  assumed  that  equilibrium 
exists  between  the  carbon  dioxide  stores  in  tissue  and  blood, 
and  that  at  equilibrium  the  concentration  in  tissue  is  one- 
half  of  the  concentration  in  blood  (4).  The  equation  for 
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CO,  which 

corresponds 

to  Eq. 

(0.5  et  ♦ 

V >r'carl 

+ q (c( 

4 

dt 

i RC02,i 


(20> 


According  to  the  simplified  reaction  scheme  presented  above, 

R„„ 


-V.i  R 


'*0,.!  R 


0. 


25\a.g 


(In  muscle) 

(In  liver)  (21) 


in  which  c ■ rate  at  which  lactic  acid  is  converted  to 
glycogen  in  the  liver 
R • respiratory  quotient  £ 0.9. 

Equations  of  change  are  required  also  for  the  lactic 
acid  produced  during  anaerobic  glycolysis.  In  this  case, 
equilibrium  does  not  exist  between  lactic  acid  in  muscle  and 
bleed.  Therefore,  we  assume  that  lactic  acid  is  transferred 
from  muscle  to  blood  at  a rate  which  is  proportional  to  the 
difference  in  concentration;  l.e., 


tb 


Equations  13  and  14  reduce  to  the  following  forms: 


4 

and 


dt 


R, 


LA,  1 


DLA(CLAt.l  ' CLAV1) 


(22) 


(23) 


sdc^ 

4 dt 


+ q(CLAb,i  _ CLAb,i-l)  " DLA(CLAt,i  “ Sj^.i) 


(24) 


The  net  rate  of  production  of  lactic  acid,  R.  , is  defined 

JR  follraj*  LA»1 


05,i 


- - 125  "lA.0.1 
in  which  c>|  - 


(in  muscle) 
(in  liver) 


(25) 


(26) 


Values  for  the  constants,  D and  k , are  not  available  in 
the  literature;  they  must  be  estimated  by 

comparing  computed  results  with  a small  number  of  reported 
experimental  observations. 

Equations  of  change  for  inert  gases  are  also  easilv  ob- 
tained. Again  it  is  assumed  that  local  equilibrium  exists 
between  gas  dissolved  in  tissue  and  blood.  Because  there 
is  no  production,  Eqs.  13  and  14  reduce  to 


£i 

b7  eft 


Wj.j)  - o (27) 

in  which  r ■ ratio  of  the  solubilities  in  tissue  and  blood. 


In  the  lungs  additional  terms  are  required  to  account  for 
exchange  of  gases  between  pulmonsry  capillaries  and  the 
alveolar  space.  NetaboLlt  reactions  are  neglected  In  the 
lunge,  and  we  have  for  both  reactive  and  inert  gases 


<r  S ♦ V 3T  + " ,ci 


ci-i»  - -Fb.,i 


(28) 


in  which  ^ ■ rate  at  which  the  component  ia  transferred 
* from  blood  to  alveolar  gas. 


Since  q ior  the  lungs  is  u relatively  large  quantity,  one  can 
neglect  accumulation  of  mass  in  the  capillaries  and  surround- 
ing tissue  and  approximate  Eq . 28  by 


q(C  - C.)  - -F. 
a 4 ba 

in  which  » concentration  in  the  pulmonary  artery, 
total  rate  of  transfer  along  a capillary.  The  method 
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In  developing  equations  of  change  for  large  arteries  and 
veins,  it  is  assumed  that  there  is  neither  reaction  nor  ex- 
change of  mass  with  surrounding  tissue.  Hence,  we  have  the 
following  equation  for  oxygen  in  the  arterial  pool  of  the 


i — element . 


dC 


02a,i 

“dt 


«* i (C0,a,»  ~ C05a,l) 


( 30) 


In  which  C , * oxygen  concentration  of  arterial  blood  in 
n.a.i  i . t h 

7 the  i — segment. 


0.a,m 


oxygen  concentration  of  blocd  entering  the 
the  arterial  pool. 


qai  • volumetric  t low  rate  of  blood  entering  tne 
arterial  pool. 

The  corresponding  equation  for  oxygen  in  the  venous  pool 


02v,i 


<vl  ^O.v.n  ” C0-v,i^ 


«»<Cn 


C.n  . )rdr 
0?v,i 


trie  fl 

the  ii—  element  from  the 
term  ^ under  the  Integral  refers  to  the 
oxvgen  2*  concentration  of  blood  flowing  into 
the  venous  pool  from  capillaries  of  the 
element . 


(31) 


in  which  q • volumetric  flow  rate  for  venous  blood  entering 
v * -5  nil)  clement,  and  the 


Corresponding  equations  for  other  components  can  be  ob- 
tained bv  merely  changing  the  identifying  subscript  in 
Eqs . 30  and  31 . 

Heat  and  Mass  Transport  in  the  Lungs 

When  evaluating  exchange  rates  in  the  respiratory  system, 
we  assume  that  the  lungs  can  be  treated  as  a we  11 -mixed 
region  of  constant  volume,  V , through  which  gas  passes  at 
a rate,  V . Since  heat  and  water  vapor  are  transferred  with 
ease  in  the  lungs.  It  is  assumed  that  expired  air  is  satu- 
rated with  water  vapor  at  the  temperature  of  the  respiratory 
tract,  T , defined  as  follows  (5): 


24.4  ♦ 0.32  T. 


(32) 


in  which  T ■ temperature  of  inspired  gas,  C.  The  rate  at 
which  sensible  heat  la  lost  through  the  respiratory  tract 
Is  easily  computed  as  follows 


V ‘ VE  CP  ,Tr  - V 

The  corresponding  rate  of  latent  lieat  loss  in 

f.H  'H  O • Vo.l 

V ’ VE  2T*  ( — 7To U27*T> 

in  which  AH  » molar  latent  heat  of  vaporization  of  water 


(31) 


(VO 


Qth»  . gases  are  exchanged  also  across  the  walla  of  the 
pulmonary  capillaries.  Sim.  tlis  process  for  a given  Ras  is 
limited  by  the  rate  at  which  gaa  la  released  In  the  <apillar- 
ien,  concentration  changes  in  the  lungs  depend  on  cardiac  out- 
put aa  well  as  the  ventilation  rate.  For  0,  the  equation  of 
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change  can  be  written  In  the  fern 


represent  outputs  from  warm  receptors  while  negative  errors 
reptesent  outputs  from  cold  receptors. 


H 


Vt  “uT*  ‘ &60F.  . ♦ Vp(P.  . - P(.  ) 

L dt  bd.O,  h 0_,i  0- 


Peripheral  afferents  from  receptors  in  the  skin  are  inte- 
grated bv  forming  weighted  sums  uf  the  signals  coming  from 
warn  and  cold  receptors. 


In  which  P, 


. and  P 

*V‘  °2 


are  the  partial  pressures  of  oxygen  in 


inspired  and  expired  air,  respectively.  The  numerical  factor, 
860,  appears  because  gas  concentrations  In  blood  are  mea- 
sured in  ml  of  dry  gas  at  STP/ml  of  blood;  a common  factor, 
B-47  (B  • pressure),  would  normally  appear  in  the  denomi- 
nator of  each  term.  Corresponding  equations  for  other  com- 
ponents car  be  written  by  merely  changing  the  subscript 
designating  the  gas. 


S.W  - VTSk.n.i  - T„«.l> 


s - r.  f4(t  - t . . .) 

ac  set,!  skln.l 


(36) 


Following  comm on  practice  we  assume  that  the  partial 
pressure  of  each  component  in  the  alveolar  gas  is  equal  to 
the  tension  of  that  gas  In  blood  us  it  leaves  the  pulmonary 
capillaries.  Strictly  speaking  equality  will  not  exist 
because  a difference  in  pressure  Is  required  to  drive  the 
exchange,  but  the  difference  is  usually  negligible.  Our 
analysis  would  nut  be  complicated  unduly  by  using  a more 
precise  expression  We  assume  also  that  2X  of  the  pulmonary 
blood  flow  is  shunted  (6). 

Evaluation  of  t »e  mass  transfer  rate  between  blood  in 
the  pulmonary  capillaries  and  alveolar  gas  requires  knowledge 
of  the  gas  tension  in  blood.  The  oxygen  dissociation 
equation  and  carbon  dioxide  buffer  equations  developed  by 
Crodlns,  Buell,  and  Bart  (7)  are  used  in  this  paper.  Both 
the  Haldane  and  the  Rohr  effects  are  included  in  this  model. 
The  bicarbonate  concentration  dec  reases  as  lactic  acid  con- 
centration increases  18).  The  discrepancy  between  the 
concentration  of  oxyhemoglobin  predicted  using  this  model  and 
generally  accepted  values  is  not  great. 


in  which  S and  S are  the  integrated  warm  and  cold  signals, 
respectively ; F aC  Is  the  weighting  factor  for  the  !.*•_ 
element;  only  positive  terms  are  Included  In  the  sum. 

Effector  signals  depend  on  the  integrated  peripheral  afferent 
signals  and  the  error  signal  generated  in  the  brain. 

Control  Functions  for  Perfusion 

One  mechanism  available  for  regulation  of  central  temper- 
ature is  control  of  the  rate  at  which  blood  perfuses  the  skin. 
As  deep  body  temperatures  rise,  blood  flow  to  the  skin  in- 
creases, thereby  increasing  skin  temperature  and  the  rate  of 
heat  transfer  to  the  environment.  When  circumstances  re- 
quire conservation  of  heat,  subcutaneous  capillarle  on- 
strict  reducing  the  rate  at  which  heat  is  brought  to  toe 
surface  by  convection.  These  responses  are  incorporated  into 
the  model  through  two  variables,  DIL  and  CON,  which  are  the 
increase  in  cardiac  output  owing  to  Increased  perfusion  of 
the  skin  and  an  index  of  the  increase  in  peripheral  circula- 
tory resistance. 


PHYSIOLOGICAL  COM  POL  EQUATIONS 


DIL  - 1.95<T.  . - T . . ) + 0.125  (S 

brain  set, brain  a 


- S > (37) 


Feedback  Control 


The  equations  of  change  presented  above  permit  one  to 
compute  temperatures  and  concentrations  when  physiological 
parameters,  such  as  local  metabolic  rates,  tissue  perfusion 
rates,  and  the  ventilation  rate,  are  known.  These  rates,  in 
turn,  depend  on  temperature,  oxygen  and  carbon  dioxide  ten- 
sions, lactic  acid  concentration,  and  level  of  exercise. 

Since  the  formulation  of  control  equations  is  still  an  a'tive 
area  of  research,  this  portion  of  the  model  is  necessarily 
somewhat  incomplete. 

Considerable  progress  has  been  made  recently  in  obtaining 
a quantitative  description  of  thermal  regulatory  mechanisms, 
especially  under  warm  conditions  (9,  10).  Using  concepts 
drawn  from  the  cheory  of  feedback  controllers,  one  separates 
the  control  system  into  three  parts.  The  first  contains 
various  thermoreceptors  which  define  the  thermal  state  of  the 
system.  The  second  consists  of  the  control  center  which 
receives  information  from  the  receptors,  integrates  it,  and 
transmits  appropriate  commands  to  the  effector  systems,  which 
constitute  the  third  part.  Effector  commands  may  be  modified 
depending  on  local  conditions  before  being  translated  into 
effector  action. 

Thermoreceptors  have  been  identified  by  either  observing 
the  effect  of  a given  temperature  change  on  quantities  such 
as  sweat  rate  and  skin  blood  flow  rate,  or  by  observing  the 
change  in  electrical  activity  in  nerve  fibers  as  the  tempera- 
ture changes.  Areas  In  which  thermosensitivity  has  been 
demonstrated  include  the  preoptic  area  and  the  anterior  hypo- 
thalamus in  the  brain,  the  abdominal  viscera,  and  the  skin. 
Other  areas  have  been  suspected,  but  since  it  is  extremely 
difficult  to  restrict  thermal  stimulation  to  a specific  area, 
the  function  of  a given  area  Is  not  easily  evaluated. 

The  output  from  each  thermoreceptor  is  generally  defined 
in  terms  of  an  "error  signal"  which  is  the  difference  between 
the  instaneous  temperature  and  the  reference  temperature,  or 
set-point,  for  a given  area.  One  can  question  the  existence 
of  a large  number  of  setpoints,  but  an  alternative  point  of 
view  1b  that  this  form  of  the  recepcor  equation  represents  a 
linearization  of  more  complex  equations  about  an  equilibrium 
state.  Some  of  the  more  recent  work  (11)  includes  a rate-ef- 
change  term  in  the  error  expression,  but  we  have  not  included 
that  In  the  present  model.  Positive  errors  are  assumed  to 


(38) 


DIL  is  measured  in  liters/mln  and  CON  is  dimensionless.  The 
blood  flow  rate  to  each  region  is  computed  by  adding  a 
fraction  of  DIL  to  the  base  flow  rate  for  the  region  and 
dividing  the  resulting  value  bv  1 + FC  CON,  in  which  FC^  is 
a constant  factor  defined  for  the  region.  It  is  assumed 
also  that  local  factors  cause  a doubling  or  halving  of  the 
perfusion  rate  for  each  6°C  change  in  the  local  temperature. 
Hence,  the  final  expression  for  the  cutaneous  blood  flow  rate 
to  the  i—  element  contains  a multiplicative  factor, 

Exp  [(T  - T e )/8.66j.  Although  the  gross  responses 

described6$Athe  S§d4l  have  been  verified  through  measure- 
ment of  the  mean  thermal  conductance  under  various  condi- 
tions, one  should  not  conclude  that  our  state  of  knowledge 
is  complete,  or  even  generally  satisfactory. 

For  example,  the  model  does  not  describe  the  observa- 
tions of  Spealman  (12)  who  measured  directly  the  blood  flow 
rate  to  the  hand  immersed  in  water  at  various  temperatures 
and  showed  that  the  perfusion  rate  declines  with  declining 
bath  temperature  to  10°C.  Below  this  temperature  there  Is  a 
striking  ini rease  in  the  hand  blood  flow  rate,  which  fre- 
quently attains  values  comparable  to  those  observed  in  a 
water  bath  at  35°C.  If  the  sublect  lb  comfortably  warm 
except  for  the  hand,  alternating  periods  of  vasodilation  and 
vasoconstriction  (hunting)  occur.  When  the  subject  is 
chilled,  the  increase  in  blood  flow  is  not  so  pronounced.  If 
the  subject  is  uncomfortably  warm,  the  blood  flow  rate  to  the 
hand  remains  high  even  in  water  at  10°C. 

The  perfusion  rate  used  in  the  model  for  muscle  depends 
on  the  local  metabolic  rate  and  the  end-capi 1 larv  oxygen 
tension.  When  Pp  45  T.  rr,  the  resting  value  is  multiplied 

bv  a factor  2 1 ♦ (45  - P )/8.  During  exercise  a 

term  which  is  proportional  to  the  difference  betvein 

the  metabolic  rates  in  the  working  and  resting  states  is 
added  to  the  resting  blood  flow  rate.  It  is  assumed  that 
perfusion  rates  do  not  change  Instantaneously,  but  instead 
change  exponentially  with  a time  constant  of  30  seconds  for 
all  tissue  except  brain  which  has  a time  constant  of  O.i 
seconds . 
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During  vigorous  exercise  blood  1 low  to  inactive  muscl<-  is 


reduced  In  order  to  nrec  the  demand  of  active  muscles  without 
exceeding  1 imitations  on  cn'dta  output  (13,  14).  This 
feature  was  lncorpor a.  ed  into  ti  e mode!  by  multiplying  each 
flow  rate  for  muscle  l>v  a ! a<  tor  l - O.uOOU  exp  (0.32  car- 
diac output);  cardfa«  output  Ij  measured  in  liters/mln. 

Since  cardiac  output  rc'resents  the  summation  of  local  per- 
fusion rates,  its  evaluation  requires  the  solution  of  a tran- 
scendental equation. 

Artec tal- venous  oxygen  difference  can  bv  used  to  ascer- 
tain the  blood  flow  rate  for  legions  which  are  supplied  by 
an  accessible  artery  and  are  drained  by  a representative 
vein.  In  this  way  it  has  been  established  that  the  cerebral 
blood  flow  rate  varies  considerably , apparently  in  response 
«.o  varying  arterial  carbon  dioxide  tension  (6).  Experimental 
data  reported  by  l.awb*rtsen  et  al  (13,  b)  were  used  to  define 
the  relationship  which  is  used  in  the  model  to  compute  the 
rate  of  cerebral  blood  flow.  Furthermore,  when  the  end- 
capillary oxygen  tension  falls  below  43  Torr,  the  base  flow 
rate  is  multiplied  by  1 + (43  - )/45.  We  assume  also  that 

high  arterial  oxygen  tensions  (P  J > 500)  reduce  the 

cerebral  blood  flow  rate  by  a factor  l - 0.00005  (P  - 500). 

aU2 

Admittedly,  the  perfusion  control  equations  Incorporated 
into  the  current  model  represent  only  a first  approximation 
to  reality.  More  precise  relationships,  such  as  those,  pre- 
sented by  Rowell  (13),  should  be  built  into  the  model  and  the. 
resulting  effect  on  system  performance  should  be  evaluated. 


Control  Functions  for  Sweat! ng 


Man's  ability  to  regulate  his  central  temperature  by 
controlling  only  cutaneous  perfusion  is  limited  to  a rather 
narrow  range  of  environmental  conditions  and  exercise  levels. 
When  heat  must  be  transferred  to  the  environment  at  higher 
rates  than  permitted  by  convection  and  radiation  alone,  the 
subject  begins  to  sweat.  Since  sweating  is  such  an  important 
thermoregulatory  mechanism  and  one  that  can  be  observed 
directly,  numerous  studies  have  been  -eported  in  the  litera- 
ture. Nevertheless,  important  questions  remain  unanswered. 


One  concerns  the  importance  ot  a non-thermal  component 
in  the  forcing  function  for  sweating  during  exercise.  Van 
Beaumont  and  Bullard  (16)  provided  evidence  in  support  of 
such  a component  with  their  observation  that  an  increase  in 
the  sweat  rate  occurs  within  a few  seconds  of  the  onset  of 
exercise  bv  a subject  who  Is  already  sweating.  The  Increase 
in  sweat  rate  occurs  prior  to  a significant  rise  in  any  body 
temperature.  Nadel  et_  al  (17)  report  that  the  rapid  increase 
in  sweating  is  transient,  persisting  only  for  the  first 
minute  or  two  of  exercise.  Thereafter,  the  rate  of  sweating 
seems  to  be  described  adequately  by  a proportional  feedback 
control  mechanism.  The  efferent  outflow,  S in  Kcai/min,  is 
computed  as  follows 


5 . 33(T.  . - T „ . , ) 

brain  sec, brain 


+ 0 . 4 8 ( S - S ) 


(39) 


A fraction,  FJt  of  the  efferent  outflow  is  added  to  the  base 
evaporation  rate  of  the  l— - element,  and  the  result  Is  multi- 
plied by  the  factor  exp  ( (T  . . . - T . 3/101  to  determine 

► u i i t I,  skin.i  set, l 

the  local  rate  of  sweating. 


The  rate  of  evaporation  may  be  smaller  than  the  rate  of 
sweat  production  if  the  environmental  humidity  is  too  high. 
Furthermore,  it  appears  that  water  standing  on  the  skin  in- 
hibits the  secretion  of  sweat  (18),  but  this  et feet  is  not 
included  in  the  model. 

Control  Functions  for  Shivering 

Although  shivering  can  be  measured  easily,  the  control 
function  for  shivering  has  remained  ill-defined  because 
steady-state  measurements  are  difficult  to  obtain.  However, 
It  seems  likely  that  both  the  skin  and  central  temperatures 
must  fall  below  their  threshold  levels  before  shivering  can 
occur.  Experimental  studies  (19)  appear  to  be  described 
adequately  by  a multiplicative  control  function  of  the  form 


S.  • 0. 35(T  . . - Tw  , )S 

h set, brain  brain  a 


(40) 


Is  measured  in  Kcal/min.  The  increase  In  metabolic  rate 
of  the  1—  element  owing  to  shivering  Is  calculated  bv 
multiplying  Sh  by  a fraction,  FCHj  . 


Con t_r ud  Equations  I • > : the  Ventilation  Rat e 

Since  breathing  frequency  and  tidal  volume  are  relatively 
easy  to  measure,  one  might  expect  that  a generally  applicable 
control  equation  would  be  available  for  ventilation.  However, 
this  is  not  the  case  because  the  primary  sensors  which  affect 
ventilation  are  inaccessible  in  man.  It  is  known  that  re- 
ceptors are  located  in  the  aortic  and  carotid  bodies,  as  well 
as  in  the  brain.  These  receptors  are  stimulated  by  changes 
in  both  pH  and  the  CO^  tension  of  arterial  blood  perfusing  the 
region  Since  changing  CO  tension  also  changes  pH,  ir  is 
difficult  to  ascertain  whether  the  receptors  are  responding  to 
the  change  in  CO.,  tension  or  hydrogen  ion  concentration.  The 
peripheral  receptors  also  respond  to  large,  changes  in  0^ 
tension;  ventilation  is  stimulated  by  hypoxia  and  depressed 
bv  0-  tensions  in  excess  of  1 atmosphere.  Another  complicat- 
ing factor  is  that  cerebral  perfusion  rate  is  affected  by 
changes  in  0-,  and  CO.,  tensions.  Furthermore,  there  exists  a 
drive  of  undetermined  origin  which  causes  a prompt  increase 
in  the  ventilation  rate  at  the  beginning  of  exercise. 

Lambert  sen  (6)  has  presented  an  excellent  discussion  of 
factors  which  must  be  considered  in  developing  a control  model 
for  ventilation.  His  view  of  ventilatory  control  is  based  In 
part  on  an  interpretation  of  step  response  data.  A thorough 
diinujsion  of  these  experiments  was  presented  recently  by 
r,t l f and  and  Lambertsen  (20),  who  concluded  that  three  sensors 
contribute  to  the  overall  response.  A fast  component,  which 
has  a delay  time  of  6-8  sec  and  a time  constant  of  5 sec  was 
assigned  to  the  peripheral  sensor.  Two  components,  which  have 
delay  i imes  of  lb-18  sec  and  time  constants  of  10  to  100  sec, 
were  assigned  to  sensors  located  in  the  central  nervous 
system. 

Further  support  for  the  "dual  center  hypothesis"  is  pro- 
vided bv  mcasui ement  of  the  steady-state  ventilatory  response 
to  a change  in  blood  pH  produced  bv  infusion  of  fixed  acids 
and  bases  while  CO^  tension  is  held  constant . These  experi- 
ments (21)  show  that  an  acid  or  base  induced  change  in  pH 
produces  approximately  30  percent  of  the  ventilatory  re- 
sponse produced  by  an  equal  CO-  induced  change  i£  pH. 

Lambert  sen  feels  that  both  centers  respond  to  |H  J,  rather 
than  P ^ « hut  one  center  lies  behind  a blood-brain  barrier 
which  2 transmits  00^  an«£  not  H , while  the  other  center  is 
responsive  to  blood  [ H ).  Although  the  magnitude  of  the 

chang|  Is  smaller,  the  ventilatory  response  to  a step  change 
in  [H  ] with  held  constant  is  just  as  rapid  as  the  re- 

sponse to  a ^2  step  change  in  p jCl)  • The  center  which 
responds  to  [H  ) must  be  vc-rv  well  ‘ 2 perfused  and  lack 
an  appreciable  diffusion  barrier  for  Ions. 

Several  attempts  have  been  made  to  stimulate  the  svstem 
by  a sinusoidally  varying  alveolar  CO,  tension.  The  most 
impressive  effort  was  reported  by  Swanson  and  bellviile  (22), 
who  recognized  that  as  the  frequency  increases,  the  relative 
importance  of  the  peripheral  receptors  increases  also.  Tn 
this  way,  they  were  able  to  differentiate  between  responses 
generated  by  the  peripheral  and  central  receptors. 

Recent  papers  by  Stoll  (23),  Bellviile,  et  al  (24),  and 
Swanson  and  Bellviile  (22)  present  independent  sets  of  fre- 
quency response  data  which  are  in  good  agreement  with  each 
other.  When  one  compute*  the  frequency  rospon-  * function 
using  the  three -component  dynamic  model  proposed  bv  fie  If  and 
and  Lambertsen,  it  is  found  that  the  predated  and  computed 
amplitude  ratios  are  ir.  reasonable  agreement,  but  t lie  pre- 
dicted phase  lag  is  mm  I)  greater  than  the  observed  lag.  The 
large  predicted  phase  lag  is  a consequence  of  the  16-18 
second  delav  tine  assigned  to  the  rwo  central  receptors;  a 
delay  time  of  6 seconds  would  have  produced  much  better  agree- 
ment. No  obvious  explanation  for  the  discrepancy  between  the 
conclusions  based  on  step  response  and  frequency  response 
data  is  available. 

Recently  Mllhorn  and  Reynolds  (25)  have  used  "exponen- 
tial peeling"  to  analyze  transient  response  data  taken  in 
their  laboratory.  Their  experience  with  this  technique,  whi.h 
is  the  one  used  by  Gelfand  and  Lambertsen,  led  them  to  ques- 
tion the  validity  of  the  control  scheme  proposed  bv 
Lambertsen  (6).  Although  It  is  impossible  at  the  present 
time  to  resolve  conclusively  the  differences  between  models 
proposed  by  various  invest igators , the  situation  is  not  hope- 
less There  is  generally  good  agreement  ibout  the  steadv- 
stat.  response  to  an  Increase  in  CO,  partial  pressure.  Since 
various  investigators  also  agree  about  the  nature  ,»f  the 
slowly  responding  sensors,  it  is  possible  to  construct  a 
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respiratory  controller  which  describes  most  rom-von  1 y o.  . urrlng 
transients  with  acceptable  accuracy.  The  controller  incc.por- 
ated  into  the  present  model  is  based  on  Lambert  sen's  work  (6), 
but  ocher  models  could  have  been  used. 


It  it»  hypothesized  that  three  hydrogen  ion  concentra- 
tions are  important.  One  is  sensed  by  a peripheral  receptor 
which  is  exposed  to  arterial  blood.  The  other  two  are  central 
receptors,  one  of  which  Is  exposed  to  blood  while  the  other 
lies  behind  a barriei  that  is  permeabh  to  carbon  dioxide  but 
not  hydrogen  ions.  We  assume  that  the  second  central  receptor 
responds  as  a first-order  system,  i.e.. 


" cs-c°. 


pcs,co,) 


(41) 


in  which  ^ - CO^  tension  at  the  central  sensor  site, 


*c°2 

k 


■ a constant  ■ 1/25  sec. 


The  three  hydrogen  ion  concentrations  art  weighted  as  follows 
co  form  a single  control  variable,  CVE. 


CVE  - 0.12  C 


a ,H'r 


°*44  Cv,H+  + 0-44  CCS.H+ 


(42) 


C - arterial  concentration  in  the  thorax  2 concentration 
* in  the  head,  and 

Cv  * end-capillary  concentration  in  the  brain. 

Tne  relationship  between  Vg  ^ and  CVE  is  shown  in  Fig.  3. 


It  is  known  that  under  certain  conditions,  the  minute 
volume  is  influenced  also  by  the  arterial  oxygen  tension.  The 
value  obtained  from  Pig.  3 is  modified  as  follows  to  form 

ve,p' 

VE  f " VE  f ' C0ER(VS.f  " 6)  (43> 

in  which  COER  * 0.i(P  - i04)/656.  If  P _ *"  37,  another 

term  is  added  to  V'  .2  a°2 

0.36  exp  (0.23(37  - Pg0  )J  (44) 

These  equations  apply  to  resting  subjects.  During  exercise  Vg 
exceeds  Vg  ^ by  a considerable  amount,  presumably  due  to  the 
influence  ’ of  neurogenic  factors  (26).  A change  in  central 
temperature  also  affects  the  ventilation  rates 


water.  When  a gas  b involved  either  internal  or  external  to 
the  garment,  one  ran  specify  whet  tier  ic  Is  helium  or  air.  In 
the  case  of  a closed  cell  foamed  neoprene  garment,  the  thermal 
conduct ivity,  density,  specific  heat,  ar.d  thickness  are 
adjusted  for  pressure. 

In  solving  the  heat  conduction  equation  numerically,  one 
ha6  a certain  amount  of  freedom  in  formulating  the  finite- 
difference  equations.  We  chose  to  use  an  implicit  set  of 
difference  equations  of  the  Crank-Nlchc isor:  type,  which  re- 
quired the  simultaneous  solution  of  sane  two  hundred  and 
fifty  equations  at  each  time  step.  Although  the  coefficient 
matrix  did  not  have  the  classical  tridiagonal  form  that  one 
customarily  obtains,  the  set  of  equations  was  still  well- 
conditioned  for. solution  by  Gaussian  elimination.  Readers 
who  are  interested  in  the  details  of  the  solution  can  find  a 
fairly  complete  discussion  in  a previously  published  paper 
(27). 

REPRESENTATIVE  RESULTS 


The  number  of  cases  that  can  be  analyzed  using  this 
program  is  nearly  infinite.  Results  for  three  typical  cases 
will  be  discussed  in  this  section. 

Exposure  to  Elevated  Temperatures  in  a Hyperbaric  Chamber 

It  is  widely  recognized  that  divers  are  subject  to 
chilling  in  a hyperbaric  chamber  containing  helium.  However, 
the  equally  great  danger  posed  bv  an  elevated  chamber  temper- 
ature was  not  fully  appreciated  until  two  divers  perished  in 
Scotland  from  accidental  hypert heimia . Therefore,  this  ib  a 
worthwhile  case  to  analyze.  Two  sets  of  computations  will  be 
discussed  - one  far  a one  hour  exposure  to  1 ata  air  at  41°C 
and  the  other  for  a similar  exposure  to  7 ata  helium  (corres- 
ponding to  a depth  of  62  meters).  In  each  case  the  diver  is 
dressed  in  a short  sleeved  shirt,  cotton  trousers,  and  shoes. 

Shown  in  Fig.  5 are  graphs  of  the  rectal  and  mean  skin 
temperatures  for  these  two  exposures.  In  air  the  diver's 
mean  skin  temperature  rises  i rora  34.2°C  to  35.4°C  during  the 
first  six  minutes  of  exposure.  Then  profuse  sweating  occurs 
and  evaporative  cooling  lowers  the  skin  temperature  to  a 
reasonable  level.  During  the  entire  exposure,  the  diver 
receives  more  heat  by  convection  than  he  loses  by  evapora- 
tion, and  the  rectal  temperature  Increases  from  37.2°C  to 
nearly  38°C.  While  this  would  be  an  unpleasant  experience, 
it  would  not  be  damaging  to  health 


SOLUTION  OF _THE  EQUATIONS 

The  equations  described  in  the  previous  section  were 
solved  numerically  using  finite-difference  techniques  This 
procedure  provides  considerable  latitude  in  dealing  with  para- 
meters and  variables,  such  as  thermal  conductivity,  specific 
heat,  metabolic  rate  and  perfusion  rate,  which  varv  with  posi- 
tion and  time.  Each  of  the  15  major  elements  is  subdivided 
into  not  more  than  8 concentric  regions,  In  each  of  which  the 
properties  are  uniform.  A representative  arrangement  Is 
shown  in  Fig.  4.  By  properly  defining  these  regions  and  as- 
signing physical  properties,  one  can  simulate  many  different 
kinds  of  thermal  systems. 

Auxilliary  programs  and  subroutines  have  been  written  to 
assist  with  the  evaluation  of  physical  properties.  Since  one 
objective  of  this  study  has  been  to  evaluate  the  accuracy  of 
the  model  by  comparing  computed  results  for  a specific  ex- 
periment with  the  corresponding  experimental  results,  it  was 
helpful  to  prepay  an  auxilliary  program  which  could  generate 
geometric  values  and  physical  properties  for  a specific  indi- 
vidual. Input  variables  for  this  program  include  the  sub- 
ject's height,  weight,  percent  fat,  and  the  length,  perimeter, 
and  coronal  and  sagittal  diameters  of  each  of  the  major  ele- 
ments. Quantities  such  as  the  percentages  of  bone,  blood,  and 
viscera  are  taken  from  published  tables.  Output  variables, 
which  are  saved  In  a file  that  can  be  read  by  the  main  pro- 
gram, include  the  locations  of  the  radial  nodes  for  each  of 
the  elements,  ldentif l cat  ion  of  the  boundary  points  which 
separate  adjacent  property  regions,  and  specification  of  mean 
properties  for  each  of  the  regions  out  to,  and  including,  the 
skin. 

A subroutine  is  available  for  evaluating  the  thermal 
properties  of  various  kinds  of  outer  garments,  including  a 
shirt,  socks,  trousers,  shorts,  dry  suit  with  woolen 
underwear,  and  a wet  suit  which  mav  be  supplied  with  warm 


The  corresponding  graphs  tor  exposure  to  helium  Indicate 
that  this  condition  is  potentially  dangerous.  Except  for  a 
short  pause  in  the  rising  skin  temperature,  both  temperatures 
Increase  steadily  during  the  entire  exposure.  Several  factors 
contribute  to  the  severity  of  the  helium  exposure.  One,  of 
course,  is  that  the  thermal  conductivity  of  helium  is  six 
times  that  of  air  which  enhances  the  rate  at  which  heat  1$ 
transferred  to  the  skin  by  convection.  Another  is  that  the 
greater  density  of  gas  that  is  breathed  promotes  heat  gain 
through  the  respiratory  tract.  The  third  factor,  which  is  not 
so  obvious.  Is  that  the  diffusivity  of  water  vapor  decreases 
with  increasing  gas  density.  This  means  that  the  maximum 
rate  of  evaporation  decreases  as  the  pressure  Increases,  and 
sweating  becomes  ineffective  as  a means  for  coping  with 
rising  temperatures. 

Evaluat ion  of  a Commonly  Used  Wet  Suit 

Wet  suits  made  from  closed  cell  foamed  neoprene  are  often 
used  to  provide  thermal  protection  during  relatively  shallow 
dives.  A typical  ensemble  consists  of  1/4  In.  thick  "Farmer 
John"  pants  and  jacket  with  hood,  a 1/8  In.  thick  vest  with 
hood.  Mu  in.  thick  soft-soled  boots,  and  three-fingered 
mittens.  A face  mask  and  fins  arc  also  worn.  The  thermal 
response  of  a diver  of  medium  build  while  wearing  this  suit 
for  two  hours  in  water  at  4,  16,  and  27°C  was  studied  at  1 
ata.  In  each  case,  the  diver  alternated  periods  of  rest  with 
periods  of  moderate  swimming  as  shown  In  Figs.  6 through  8. 

It  appears  that  this  suit  Is  quite  adequate  for  use  in 
16°C  water.  There  Is  a modest  decrease  In  rectal  temperature 
during  the  resting  period,  but  the  trend  is  reversed  during 
exercise.  It  is  interesting  to  note  that  there  Is  an  under- 
shoot following  the  end  of  the  rest  period  and  an  overshoot 
following  the  exercise  period.  The  rate  of  gas  consumption 
during  the  exercise  period*  is  approximately  double  the 


48* 


resting  rale.  As  the  dive  prognoses,  there  i -•  .»  long  term 
trend  toward  in.  rt-asiag  gas  consumption,  which  can  be  attri- 
buted to  low  level  shivering  cauoe  1 by  the  decrease  In  rectal 
and  skin  temperatures.  Nevertheless,  it  appears  that  the 
diver  is  quite  rafurtable  under  these  conditions. 

Corresponding  graphs  obtained  when  the  water  temperature 
is  4°C  indicate  that  the  s lit  does  not  provide  adequate  pro- 
tection for  long  periods  of  exposure  at  this  temperature. 

The  rectal  temperature  falls  1°C  during  the  first  30  minutes 
>f  the  dive.  This  results  in  vigorous  shivering  with  a 
corresponding  increase  in  the  race  of  gas  consumption.  Some 
divers,  especially  those  who  ate  some whet  falter  than  the 
model  diver,  could  tolerate  this  exposure  for  an  extended 
period  of  time,  but  they  would  be  quite  uncomfortable  and 
fatigued  at  the  conclusion  of  the  dive. 

A third  set  of  computations  was  performed  for  2.7°C  water. 
The  results  are  shown  in  Fig.  3,  where  one  car.  see  that  the 
diver  overheats  during  the  exercise  period.  From  these 
results  one  can  conclude  that  this  suit  is  ideal  for  shallow 
dives  when  the  water  temperature  is  approximately  16  C. 

The  effect  of  depth  on  performance  while  wearing  this 
suit  is  illustrated  in  Fig.  9 where  graphs  are  presented  for 
a 30  meter  dive  in  16°C  water.  It  is  obvious  that  the  suit 
is  very  inadequate  under  these  conditions.  Violent  shiver- 
ing occurs  early  in  the  dive,  and  it  is  doubtful  whether  the 
diver  could  remain  in  the  water  very  long.  The  flat  gas 
consumpr  ion  rate  during  exercise  occurred  because  of  an 
arbitrary  maximum  value  specified  in  the  program. 

The  1 at.*  results  are  consistent  with  observations  made 
at  the  Naval  Clascal  Systems  Laboratory  in  Panama  Lily, 

Florida  (23).  When  the  wet  suit  used  in  these  calculations 
was  tested,  it  was  found  chat  the  suit  was  ideal  for  use  in 
16  (.*  water.  At  4 C shivering  occurred  and  many  divers  were 
unable  to  complete  a 6 hour  dive.  The  2 7°C  dives  were  con- 
ducted wearing  a partial  suit  to  avoid  overheating.  Hence, 
tne  model  does  appear  to  be  capable  of  simulating  reasonably 
well  the  performance  of  divers  under  these  condition?.. 

Use  of  a Liquid  Heated  Dry  Suit 

Prolonged  dives  depths  below  200  meters  require  the 
use  of  auxiliary  heating  to  keep  the  diver  comfortable.  Al- 
though they  have  not  been  used  to  date,  closed-circuit  liquid 
heated  suits  would  appear  to  be  more  efficient  than  the  open- 
circuit  suits  that  are  generally  used.  Therefore,  several 
cases  were  analyzed  to  determine  whether  the  closed-circuit 
concept  Is  reasonable. 
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1.  All  Q's  are  measured  in  Real /min. 

2.  All  T's  are  measured  in  °C. 

3.  W is  measured  in  kgm/min. 

4 . Notat ion  - 

metabolic  rate  not  including  shivering 


MET 

Wuq 


CA-S 

Qppj. 


flow  rate  through  liquid  heated  garment 
inlet  temperature  to  liquid  heated  garment 

temperature  rise  across  regenerative  heat 
exchanger 
rectal  temperature 

temperature  of  brain 
temperature  of  chest  skin 
temperature  oi  head  skin 
temperature  of  calf  skin 

rate  ot  heat  loss  through  the  respiratory  tract 
rate  of  heat  loss  by  convection 

rate  of  heat  generation  due  to  shivering  (nega- 
tive values  denote  sweat  production,  without 
rvaporat ion) 

rate  of  heat  input  bv  liquid  heated  garment 


In  the  present  version  of  the  model,  a liquid  heated  gar- 
ment can  be  defined  by  specifying  the  heat  transfer  coeffi- 
cient and  area  of  contact  for  each  element  of  a suit  such  as 
the  one  illustrated  in  Fig.  10.  Since  the  fluid  which  serves 
as  the  heat  source  may  flow  serially  past  several  heat  ex- 
change areas,  it  is  necessary  to  specify  the  flow  path  for 
each  stream.  This  is  done  bv  indexing  each  segment  of  the 
path  and  developing  a table  In  which  are  specified  (1)  the 
anatomical  element  with  which  each  exchanger  segment  Is  in 
contact,  (2)  the  indexes  of  adjacent  segments  which  feed  into 
each  segment,  and  (3)  the  fraction  of  the  total  flow  which 
passes  through  each  segment.  The  structure  of  the  program 
allows  one  to  alter  a liquid  heated  garment  by  changing  a few 
data  statements.  Hence,  a situation  such  as  loss  of  flow  to 
a particular  segment  owing  to  obstruction  of  the  tubing  can 
be  simulated  with  ease. 

To  Illustrate  the  use  of  a liquid  heated  suit,  calcula- 
tions were  performed  for  a diver  wearing  a 0.3  in.  thick 
foamed  neoprene  dry  suit  over  a liquid  heated  garment  having 
a thickness  of  1/8  inch.  In  some  cases  the  breathing  gas 
passes  through  a regenerative  heat  exchanger  which  produces  i 
25°C  difference  between  the  temperature  of  the  supply  tank 
and  the  temperature  at  the  mouth.  Humidification  of  the 
inspired  gas  also  occurs  in  the  heat  exchanger.  When  the  re- 
generator Is  not  in  the  circuit,  the  temperature  and  dew 
point  of  the  inspired  gas  are  set  equal  to  the  water  tempera- 
ture. Temperatures  and  heat  transfer  rates  were  computed 
for  a 50  minute  exposure  to  0°C  water  at  a pressure  of 
10  ata.  Computed  results  are  summarized  in  Table  1. 


The  results  Illustrate  several  points.  One  Is  that  this 
environment  presents  a severe  thermal  stress  for  the  diver. 

The  rate  of  heat  transfer  to  the  environment  depends  on  the 
rate  at  which  warm  water  circulates  through  the  liquid  heated 
garment,  but  it  is  never  less  than  7 Kcal/min  (or  the  cases 
that  were  studied.  Therefore,  auxiliary  heating  is  required, 
except  when  the  diver  is  working  vigorously  (metabolic  rate  * 

6 Kcal/min).  Even  for  the  working  diver,  heat  loss  through 
the  respiratory  svstem  must  be  limited  if  he  is  to  maintain 
his  thermal  balance. 

Calculations  performed  assuming  that  the  breathing  cir- 
cuit contains  a regenerative  heat  exchanger  show  that  this  is 
a very  valuable  device.  It  is  capable  oi  reducing  the  rate  of 
heat  loss  through  the  respiratory  system  for  the  working  diver 
from  2.6  Kcal/min  to  0.8  Kcal/min.  This  is  especially  valu- 
able because  heat  lost  through  the  respiratory  tract  comes 
from  deeper  tissues  which  are  well  perfused  with  blood.  Fui- 
thermore,  the  device  is  passive  and  requires  no  heat  input, 
although  it  might  be  worthwhile  to  add  heat  if  the  device 
were  to  be  used  at  greater  depths. 

When  heat  loss  through  the  respiratory  system  is  not  re- 
stricted, it  appears  that  a considerable  difference  may  deve- 
lop between  rectal  and  brain  temperatures.  Indeed,  the  rectal 
temperature  may  increase  slightly  owing  to  heating  by  the 
venous  return  from  working  leg  muscles  while  the  temperature 
of  the  brain  decreases  owing  to  loss  of  heat  through  the 
respiratory  tract.  Experimental  confirmation  of  this  possi- 
bility would  lend  credence  to  the  model. 
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Finally,  it  should  bo  noted  that  *ome  tom  ol  thermal 
control  will  be  required  for  the  liquid  heated  garment . i>ne 
can  vary  either  the  flow  rate  or  temperature  ot  the  Incoming 


hot  water,  but  only  variation  of  t law  rate  was  Invest igated  8. 

in  this  aeries  of  >:a]  r uiatlons . The  inlet  temperature  wans 
held  <011313111  «;  *4.4°C. 

Whet  the  diver  works  at  .1  Metabolic  rate  of  6 Kcai/mln,  9. 

he  will  be  comfortable  with  very  tittle  auxiliary  heating. 

Indeed,  a water  flow  rate  of  C.0’6  Kgm/min  produces  some 
sweating,  which  is  uodesi 1 abl e . However,  the  'eating  diver 
will  cool  rapldlv  u:.der  the  same,  conditions.  When  the  meia-  10. 

bolic  rate  is  reduced  to  !.82  Real /min,  centra,  temperatures 
begin  to  decline;  and  a.  ter  10  minut-s,  the  diver  is  shiver- 
ing. Fifty  minutes  of  exposure  results  in  a rail  in 

brain  temperature , whi.-h  causes  severe  shivering.  If  the  11. 

flow  -ite  of  hot  water  is  Increase,  to  0.38  kgm.ir.in,  b.21 
K al/min  are  released  by  re  liquid  heated  garment-  This  is 
4./  KcaWttln  greater  than  the  rate  of  heat  transfer  at  the 
lower  flow  rate.  However,  the  rate  of  heat  loss  to  the  12. 

environment  also  increases  from  7.5  to  9.1  Kcal/mln.  Hence, 

■ Ire  thermal  1 id  on  the  diver  Is  redu> . d by  approximately 
1 Kcal/mln.  His  central  cempeiature  still  falls  and  he  still 
shivers,  but  not  as  severely  <ij  at  the  1 ?w  flow  rate.  13. 

One  •oul*'  conceivably  develop  a passive  suit  that  pro- 
vide adequate  thermal  protection  for  the  resting  diver,  but 
su»  h .*  suit  would  . ..use  the  working  diver  to  sweat.  Vasa-  14. 


motor  control  alone  dees  not  allow  man  t«  adapt  to  a wide 
range  of  cond.tlons,  and  ocher  mechanisms  must  be  invoked. 

For  the  diver  tnls  could  be  -ontrol  of  the  flow  rate  to  the 
liquid  heated  garment 

15. 

CONCLUSIONS 

Even  tc.iugn  there  are  ...ops  in  our  knowledge  of  the  human 
therr.il  system.  It  is  possible  co  construct  a mathematical 
model  wrsi<  h destribes  many  Important  features  of  thermal  re- 
sponse In  man.  By  combining  a reasonable  physiological  model 
with  invent!  al  engineer  in;,  analysis,  it  appears  that  one 
can  simulate  the  thermal  behavior  of  divers  under  a variety 
of  conditions.  The  next  step  Is  to  use  the  model  under  care- 
fully controlled  conditions  s->  that  shortcomings  can  be 
corrected  ar.-l  confluence  developed  in -the  model.  If  this  can 
be  accomplished,  cur  ability  to  design  improved  diving 
systems  should  bo  greatly  enhanced. 
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Schematic  diagram  showing  the 
geometric  arrangement  of  major 
elements  and  the  circulatory 
svstem. 


Hemoglobin,  c 


Myoglobin,  c. 


Fig.  2.  Dissociation  curves  for  myoglobin  and  hemoglobin 


Relationship  between  the  ventilation  rate  and  the 
integrated  signal  from  ohemoreceptors . 


1 . LARGE  ARTERIES  AND  VEINS 

0 TISSUE  , SMALLER  ARTERIES  AND  VEINS 
BT  . MUSCLE 

IB  SUBCUTANEOUS  FAT  AND  SKIN 

2 LIQUID  HEATED  GARMENT 
H 6AS  SR&CE 

SII  INSULATING  OUTER  GARMENT 


Typical  arrangement  of  radial 
element  such  as  the  forearm. 
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Computed  rectal  temperature  and  minute 
volume  during  a shallow  dive  in  water  at 
4°C  with  two  30  minute  exercise  periods. 


Computed  rectal  and  mear.  skin  temperatures 
during  exposure  to  1 ata  air  and  7 ata 
helium  at  41°C. 
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Computed  rectal  temperature  and  minute 
volume  during  a shallow  dive  in  water 
at  16°C  with  two  30  minute  exercise  periods 


Computed  rectal  temperature  and  minute 
volume  during  a shallow  dive  in  water  at 
2 7 °C  with  one  30  minute  exercise  period. 
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